An energetic model for macromolecules unraveling 
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We propose a simple approach, based on the minimization of the total (entropic plus unfolding) 
energy of a two-state system, describing the stretch-induced unfolding of macromolecules (proteins, 
silks, nanopolymers, DNA/RNA). The model is fully analytical and enlightens the role of the differ- 
ent energetic components regulating the unfolding evolution. As an explicit example of application 
we compare the analytical results with the titin Atomic Force Microscopy experiments showing 
the ability of the model to quantitatively reproduce the mechanical behavior of macromolecules 
unfolding. 
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I. INTRODUCTION 

The last decade has shown a significant theoretical 
and experimental effort in the analysis on the thermo- 
mechanical behavior of macromolecular materials such 
as muscle tissues, spider silks [13 [4J], polymers and 
biop olymers |22l |. polysoaps [j] (5| and silks in general 
[5q . A common property (see e.g. [23, [2a], [23], |13 ]. 
and [53]) is that their macroscopic history dependent, 
dissipative properties descend by complex "semicrys- 
taUine"microstructures, constituted by flexible (poly- 
meric or protein) macromolecules reinforced by strong 
and stiff crystals e.g. in form of fillers or /3-sheets (see 
for example [1 [H). 

At small strain the stiffness of these material is mainly 
regulated by the hard fraction (secondary structure) . Un- 
der stretching a reversible hard-soft transition (e.g. /3- 
sheets unfolding in protein or crosslinks breakage in poly- 
mers) is observed. This transition has the two important 
effects of dissipating energy , due to the transition itself, 
and of a variation of the microstructure, leading to vari- 
able contour lengths of the chains. At larger stretches, 
the behavior is regulated by the entropic hardening of the 
macromolecules (primary structure). As a result, macro- 
scopically these materials shows a hysteretic behavior, 
reminiscent of the pseudoelastic behaviour characteriz- 
ing materials such as Shape Memory Alloys (see [4^] and 
references). 

The deduction of predictive models that connect the 
mesoscale properties with the macroscopic material re- 
sponse are crucial not only to describe the behavior of 
such important materials, but also in the perspective 
of the design of new bioinspircd or reconstructed bio- 
logical materials. As a consequence, an intense exper- 
imental, numerical, and theoretical effort has been re- 
cently devoted in this field [8J. From an experimen- 
tal point of view, a great impulse in this direction has 
been delivered by new experimental techniques |54{ , such 
as Atomic Force Microscopy (AFM) [S^], laser optical 
tweezers [27| . magnetic tweezers and single molecule flu- 



orescence techniques. The typical experiments is a me- 
chanically induced unfolding of a macromolecule com- 
posed of n unraveling domains, such as a polymeric 
polypeptide, dextan [50], silks 56], proteins (see [54] for 
the titin), DNA/RNA strands [53]. 

From a theoretical point of view the thermo- 
mcchanical behavior of multidomains proteins has been 
undertaken following different approaches: Molecular 
Dynamics, off lattice models, all atom Montecarlo ap- 
proaches (see [3^ and references therein), phenomeno- 
logical approaches, Statistical Mechanics energy land- 
scape analyses, with funneling [24| and Inherent Struc- 
ture models [40] ■ Molecular Dynamics theories [33] have 
been restricted by the computational effort required to 
describe the unfolding of such large macromolecules at 
the AFM loading time scale. On the other hand, the sta- 
tistical approaches for the discrete chain have been essen- 
tially based on numerical techniques, whereas analytical 
results have been obtained only in the thermodynamic 
limit hypothesis that hides the crucial role of finite size 
and discreteness of the unfolding phenomenon. 

To fix the ideas, in Fig|l]we schematically show a typ- 
ical AFM single molecule stretching experiment on an 
engineering reconstructed macromolecule domain. The 
figure shows the typical sawtooth force-length diagram, 
that can be explained as a stick-slip dynamical evolution 
in a wiggly energy landscape, characterized by multiple 
energy wells, each corresponding to a given hard-soft mi- 
crostructure configuration. Thus, for growing assigned 
end-to-end length (see [20]) the chains alternates 'slow' 
(intrabasin) steps in energy wells at fixed folded/unfolded 
configuration, followed by 'fast' (interbasins) transitions 
corresponding to the unfolding of (tipically single) crys- 
tals. These transitions are signaled by the periodic local- 
ized force drops induced by the entropy jumps due to the 
creation of new free monomers at each /3-sheet unfolding. 

The Molecular Dynamics approach in [29[, based on 
the Inherent Structure formalism, clarifies that the ob- 
served unfolding is regulated by three main timcscales: 
loading time Tioad, intrabasin relaxation time Tintra, and 
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FIG. 1: Energetic interpretation of a typical unfolding ex- 
perimental stick-slip AFM force-elongation curve (continuous 
line) and scheme of the energetic decoupling of the external 
work into unfolding (dissipated) energy Qi,i = 1,2,3,4, and 
elastic stored energy $e (dashed region). Dashed lines rep- 
resent approximating WLC curves each characterized by a 
different contour length Lc{i), i — 1,...,4; Ic is the (fixed) 
contour length increase at each unfolding event. 



interbasins transition time Tinter- Our theoretical model 
is based on the time scale separation hypothesis 
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(1) 



In this time scale regime [20| unfolding results as an alter- 
nated sequence of purely elastic, intrabasin, stick evolu- 
tions and purely dissipative, interbasins, slip transitions, 
localized at fixed unfolding length thresholds. 

In the case of AFM induced unfolding, the non dis- 
sipative hypothesis of the intrabasin evolution is sup- 
ported by the observation [13, |23 that during the (slip) 
evolution at fixed folded/unfolded configuration, the be- 
havior is fully reversible, thus indicating that the system 
relaxes to the local energy minimum. On the other hand, 
the fully dissipative hypothesis of the 'fast' (stick) 
interbasins transitions, results by the observations (see 
FigHJ that at the AFM loading time scale the /3-sheet 
unfolding events are (mainly) localized at fixed macro- 
molecule end-to-end lengths. 

We remark that the time scale separation ^ has been 
successfully adopted in other stick-slip evolutions associ- 
ated to abrupt microstructure transitions in multyvalley 
energy landscapes at 'low loading rates' and 'low temper- 
ature regimes' such as dcpinning or nucleation of new de- 
fects, dislocations and Frank-Read sources in metal plas- 
ticity (see [4g| and references therein), and Barkhauscn 
jumps in ferromagnetism [3|. 

Based on energy conservation (Gibbs equation), un- 
der our hypothesis ([1]), the external work (see [4^ for 
a theoretical discussion) can be decomposed as follows. 
By focusing again on the AFM experiment in Fig|l] sup- 
pose that we begin stretching the macromoleculc from its 
natural state. The macromoleculc follows elastically the 
first equilibrium path (0-A in FigllJ with the external 
work W accumulated as elastic energy $e (AVF = A<I>e)- 
At the first /3-sheet unfolding (path OA) - here approxi- 
mated as instantaneous with no external work W - there 



is an internal energy discontinuity []<I>[]j^ (area 0-A-B) 
that by energy conservation equals the unfolding energy 
Qi of the first /3 sheet. Similar considerations can be ex- 
tended to the next elastic and dissipative steps, so that 
the total dissipation is Q = ^ Qj = S l*^li- 

Based on previous considerations, here we consider an 
energetic approach for a two-states material, inspired by 
the model in [ll|, ll2[, where the authors describe the 
hysteresis of filled polymers and spider silks. Similarly, 
a two-state energetic approach was proposed in [9| to 
describe the helix— >-coil transition of polypeptide chains 
regulating the damage of multi-block copolymers. In this 
work the authors obtained an analytic solution in the 
thermodynamic limit of large number of breakable links. 

In the field of protein mechanics, based on the ap- 
proach in [9|, an important step in the comprehension 
of the energy competition between the unfolding and cn- 
tropic energy terms, has been delivered in j33| . In this 
paper, the author models the unfolding of a biomacro- 
molecule as a chain composed of folded and unfolded do- 
mains, both elastic with Gaussian type response, com- 
bined with an Ising-like unfolding energy. The resulting 
MD simulation well describes the unfolding effect in a 
protein macromoleculc, whereas analytical results are ob- 
tained only in the thermodynamic limit of many folded 
domains [33| . 

More recently, in [5^ a Statistical Mechanics based 
model for the stretching of titin proteins has been pro- 
posed, considering also the influence of the AFM loading 
device. A simple Ising model, neglecting the elasticity 
of both folded and unfolded domains in protein macro- 
molecules, was instead proposed in [2^ to describe the 
statistics of unfolding events. The described competition 
between the entropic energy of the unfolded fraction and 
the /3-sheets unfolding energy has been analyzed in |5l| 
and [5a] via Monte Carlo simulations combining a Worm 
Like Chain (WLC) with a two-state Bell type model for 
the unfolding. Finally, we recall the fully phenomenolog- 
ical continuum approaches recently proposed in (2| and 
[49J where the authors show the possibility of describ- 
ing the protein unfolding as an energy minimization of a 
continuum system with a non-convex internal energy. In 
particular in [2|, based on the general approach of [4y| 
and [43 for the description of the mechanical behavior of 
a bi-stable discrete chains, the authors obtain an inter- 
esting characterization of an optimality condition of the 
number of /3-shcet domains with respect to the toughness 
of the macromoleculc. 

Here, based on previous hypotheses, we obtain a fully 
analytic ratc-indcpendcnt dissipative lattice, describing 
the stretch-induced unfolding of macromoleculcs. Inter- 
estingly we show the existence of a fundamental exper- 
imentally measurable non-dimensional parameter ^ de- 
fined in p^ . representing the ratio of the elastic and 
unfolding energy of the single folded domains, that regu- 
lates the dissipation and the unfolding thresholds of each 
folded/unfolded phase configuration. The model is also 
extended to describe the possibility of the existence of 



a hierarchy of variable unfolding energies of the differ- 
ent folded domains. As we show, the effective distribu- 
tion of unfolding energies, can be deduced on the base of 
the previously described energetic analysis, using the ex- 
perimental force-elongation diagrams. Interestingly, the 
analysis of the behavior of different macromolecules sug- 
gest simple phenomcnological linear distributions. 

Aimed at the possible deduction of a three-dimensional 
continuum extension of the proposed model for biologi- 
cal tissues, that will be the subject of our future work, 
we also deduce the continuum limit of the proposed dis- 
crete model. Indeed, following statistical approaches, 
as proposed in [48|, the continuum limit model consti- 
tutes the base of a multiscale model for the description of 
the mechanical behavior of networks of modular macro- 
molecules. The main advantage of our theoretical ap- 
proach is that the behavior of the single chain also in the 
continuum limit is regulated by the experimentally mea- 
surable parameter ^ that fully characterizes the unfolding 
behavior of the macromolecule. 

Finally, as an explicit example, we focus on the AFM 
experiments of titin unfolding. As we show, under a sim- 
plifying assumption of rate-dependent effective unfolding 
energies and on the base of the phenomcnological law of 
the unfolding energy hierarchy recalled above, the analyt- 
ical model quantitatively well describes the experimental 
behavior. 



II. ENERGETIC ASSUMPTIONS 

Because in macromolecules /3-shects crystals the un- 
folding is typically an all-or-none transition, as confirmed 
also from the size of periodicity of the experimental un- 
folding lengths [52[, we model the molecule as a discrete 
lattice of n two-states (rigid-folded/entropic-unfolded) 
links (see the scheme in FigllJ. The folded/unfolded state 
of the chain is assigned by a set of internal variables Xi j 
i = 1, ...,n, such that Xi =0 {xi = 1) denotes the folded 
(unfolded) state. Thus in particular n„ = X]i=i Xi is the 
number of unfolded elements and Uf ~ n — n^ is the 
number of folded elements. 

As in the case of Fr eely Jointed Chain or Worm Like 
Chain models (see e.g. [lO[), we characterize the behavior 
of each unfolded link through its contour length Ic and 
end-to-end length I, with a free energy density (energy 
per unit length) ipf, = ipf,{ri), where 7] :~ j^ represents a 
strain measure. We assume then the limit extensibility 
condition (ySg —> oo as 77 — > 1. 

By neglecting non-local interactions (weak inter- 
action hypothesis), the total elastic energy $e ~ 

Sr=i Xi^i^e ( jt ) can be simply expressed as 



Here, 



$e = Lcipeifl)- 



V = viL,nu) 



(2) 



is the strain in the unfolded domain, 

i 

is (by neglecting the extension of the unfolded domains) 
the total end-to-end length, and 



Lc = Lc{nu) = io + nJc 



(3) 



is the total contour length of the chain. In ([3]) io de- 
notes the 'initial' (virgin) contour length of the unfolded 
domain and lc the contour length of each unfolded do- 
main. 

Indeed at equilibrium, by neglecting the elastic energy 
of the (rigid) folded domains, the total elastic energy is 
given by 4>e = Yn^iXdWei^-fi-), where l^ and k are the 
(possibly variable) contour length and end-to-end length 
of the i-th unfolded link, i.e. for all i = l,...,n with 
Xi = 1. Under an equilibrium hypothesis we have a con- 
stant stress for all unfolded links, i.e. l^-jj-^^^^ = ^■ 
Thus, for a convex energy density (monotonic derivative 
dife/dl), such as WLC or FLC, the strain is homogeneous 
in all unfolded elements, i.e. 77^ = jf = 77 = -^, for all 

i = 1, ..., ri with Xi = 1. Thus we have <I>'^ = X) ^I'P'^iv) = 

Following [^ and [33 here we consider an Ising type 
unfolding energy 



$* 



-EtiiQ-J)i^ 
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JJ2^=l (1 -X»)(l -Xi+i) ^Q{n-nu) + Jubj, 



Lc{nu) 



depending on the internal variables Xi ^-nd the number 
ribj of contiguous folded blocks in the folded/unfolded 
configuration. Here Q is the unfolding energy for a single 
domain and J is a penalizing 'interfacial' energy term 
(measuring the loss of internal energy due to the unbind 
terminal i/-bonds of each contiguous folded domain [g] ) . 
To get the total energy ^tot = -fcBT'ln[p(L,n„,7if,^.)] 
(where T is the temperature and ks is the Boltzmann 
constant), we have to know the probability p(L,ri„, n;,, ) 
of a given elastic configuration of the chain with a mi- 
crostructure corresponding to n^ and njy,. In particular, 
we have p{L,nu,nbj) = Vt{nu,nbf)pe[L,nu)ptr{nu,nbj), 
where f2(7T,„,nf,j.) represents the number of sequences 

with assigned 7t,„ and tt,;,^. , Pe[L,n^) ^ exp( °^ 'j^"'^ ) 

represents the probability of attaining a length L at given 
riu and Ptrinu,nbj:) ^ cxp( fj^) is the probabil- 
ity of a state with assigned n„. So, we obtain ^tot = 
^e{L,nu) + <&*,-("«) - TS{nu,nbj) where S{nu,nb^) = 
ks lnf7(n„, nb.) represents the mixing entropy term. 

Observe that the coupling energy term J penalizes the 
multiplicity of folded blocks, whereas the mixing entropy 
term induces multi domains configurations. In the fol- 
lowing we assume, as in [9[, that the penalizing term J 
dominates this effect, so that we always consider single 
folded domains configurations {i.e. we assume rif,, = 1, 



known as di-block approximation). This hypothesis is 
supported by the MD simulations |23| showing an un- 
folding strategy with always one single connected internal 
unfolded domain inside two boundary folded domains. 

Under these hypotheses we obtain the simple expres- 
sion of the total energy 



10^ 



$. 



$e(L, n^) + n^Q + const. 



(4) 



We remark that to avoid the introduced di-block ap- 
proximation, not always experimentally verified, one 
needs to evaluate the partition function {e.g. [33, |58| ) 
and only numerical results in the discrete model can be 
obtained. Moreover stochastic processes considering fluc- 
tuations in both the unfolding forces [35| and the un- 
folding lengths, are possible extensions of the proposed 
model. Also these extensions require the employment of 
numerical approaches. 



III. ENERGY MINIMIZATION 



Consider a WLC force-length relation proposed in (37 



fcsT 4Lp ^(1-77)2+'' 



corresponding to an energy density 



Veiv) 



ksT 
4L„ 



1 — 77 



2r,' 



(5) 



(6) 



Thus, using ([2|), the adimensionalizcd total entropic en- 
ergy of the unfolded fraction can be written as 

$e := ^ = — (^^ + 27^A LJuu). (7) 
keT ALp\l-i] 'J ^ ' ^ ' 

In order to attain analytical solutions, we consider the 
following simplified expression of the WLC energy den- 
sity: 



Veir]) 



knT 



7]^ 



4Lp V 1 - V 



(8) 



Observe that this approximation keeps the same asymp- 
totic behavior as I ^ Ic of the WLC model in ([7]). FigE] 
shows (in a log scale, stressing the differences at low val- 
ues of the force) that, while for low forces {F < 10~1 
pN) the introduced approximation is significant (as com- 
pared with the approximation in [37[), for larger forces 
the approximation is of the same order of |37| . Since in 
the low force regime the elasticity is mainly regulated by 
the PEVK and tertiary structure elasticity (see [GjI and 
[23| for details), this approximation appears inessential 
in both the qualitative and the quantitative analysis of 
the behavior during the large-forces unfolding regiine of 
interest for titin unfolding. Moreover, we remark that the 
approximation ([S|) has been shown to be inefficient in the 
low force regime in |15l |. where the authors introduce a 
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FIG. 2: Force-strain curve (log scale stresses the differences in 
the low force regime) for the WLC compared with the usual 
Marko and Siggia approximation in [33 and with the simpli- 
fied model in for Lp = 0.42. Observe that this approxi- 
mation keeps the same asymptotic behavior as / — >■ ic of the 
WLC model in 0. For low forces [F < 10^1 pN) the in- 
troduced approximation is significant as compared with the 
approximation in [37l|; f or larger forces the approximation is 
of the same order of [33 . 



Mooney Rivlin type correction to the WLC constitutive 
law. 

Thus, using ((2|) and ((3|) the total elastic energy is 

4Lp \~ri 

and, correspondingly, the total force-deformation relation 
is 



F{ri,n^ 



ksT 2i] - rf' 



4ip (1 - 77)2 ■ 
Finally, according with ((4|), the total energy is 
keT rf 



$. 



4ip 1 - 77 



Lciriu) + Qnu. 



(9) 



(10) 
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FIG. 3: Scheme of the energy minimization: with bold line 
we represent stable (global energy minimum) solutions. 

We follow a Griffith- like approach [2l[ , minimizing the 
total unfolding (fracture) energy plus elastic (entropic) 



energy, and based on ([T]) we assume that the observed 
solutions are the global minima of ^tot in (HH])- 

The first important step is to justify the experimental 
observation that the /3-sheets unfold one at a time, result- 
ing in a constant increase of the contour length (see j53|'). 
To obtain this result, we begin by evaluating the solution 
of the equations $tot(i,ri„ + rn) — <^tot{L,nu) = and 
get the intersection lengths L = L{nu,'m) (see Figl3|). 
The searched result follows by the observation that 

— g^'™ > so that if n„ is the branch corresponding to 
the global minimum, by increasing L it looses its global 
stability at the intersection with the equilibrium branch 
fiu + i (see Figl3]). Thus the chain unfolds with a sequence 
of single /3-sheets unfolding at the threshold assigned by 

$e(?7, Tiu) - $e(?7, riu + 1) = Q, 71„ = 0, ..., H - 1,: 



Lu(nu) = L(n„, 1) = _ ^ . 

(11) 
(in this formula and in the following we omit the n„ de- 
pendence of Lc). 

In (jlip we introduced the main non-dimensional pa- 
rameter of the model £, 



c 



8Lp Q 



(12) 



representing a measure of the ratio between the elastic 
and fracture energy of the single /3-sheet. Indeed we 
observe that according with (j6]) we have that "g^ " = 

ipe{^)lc measuring the elastic energy of a single domain 
when the deformation is a half of the maximum elonga- 
tion (contour length). 

It is easy to verify that L„ € (0,Lc) and that ^^ > 
0. As a result, the n„ branch corresponds to the global 
energy minimum for 

L e {Lu{nu - 1), Lu{nu)), Uu G (1, n - 1), 

representing the existence domain of the n„ branch under 
our energy minimization hypothesis. In the special cases 
of the virgin curve, with n„ = 0, we have L S (0, -L„(0)) 
and of the fully unfolded chain, with Uu ~ n, we have 
L E {Lu{n — l),ir)i where L^ is the fracture threshold 
of the fully unfolded chain. 

Using ([9]), we get the unfolding force Fu = 
F{L.u/Lc,nu) 



Fu= 



knTi 



2(e-l)ic 



4ip \ Uc i + 2Lc~ ^{Ici + ^Lcf+Ai^-l)^ 




Observe that using PT|) and ((T^ it is also possible to 
obtain an explicit relation between the unfolding forces 
and the unfolding end-to-end lengths 



Fu{Lu) 



kf 



tL.Jl.u + ./8liK+^c-Ic 



{V8LlKTl!'l, 



(14) 



It is important to observe that since ^^ < 0, under the 
hypothesis of fixed unfolding energy Q of the different /3 
sheets, the system shows a softening behavior during the 
unfolding. 
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FIG. 4: Unfolding behavior for a system of n = 6 initial 
folded domains. Here we considered the parameters: lo — 58 
nm, h = 28.43 nm, Ip = 0.36 nm, Q = 770, AQ = 420. Each 
equilibrium path is labelled by the number Uu of unfolded 
domains. 



The stretch induced unfolding of the system is shown 
with bold line in FigH) As described before, the system 
reproduces the typical experimental behavior of unravel- 
ling macromoleculcs with a regularly spaced sequence of 
unfolding events of the hard domains. Thus if we start 
loading from the virgin configuration (n„ = 0, point O in 
the figure) , the system follows elastically the equilibrium 
curve n„ = 0; until the unfolding energy Q equals the 
jump of the entropic energy due to the transition from 
the branch n„ = to the branch riu — 1 (path A-a in 
the figure) . By increasing further the assigned length the 
system follows the new branch until another sudden tran- 
sition to the branch ?t,„ = 2 is observed when it becomes 
energetically favorable (path B-b). Similar transitions 
with single domains unfoldings are then observed, until 
all the crystals unfold and the system shows a harden- 
ing behavior due to the entropic elasticity of the fully 
unfolded chain (curve f-G in the figure). 

In particular we observe that the theoretical model 
shows a softening behavior during the unfolding regime, 
with the unfolding force decreasing with the number n„ 
of unfolded domains. We remark that the experimental 
behavior of stretch induced unfolding of macromolecules 
show a variable behavior, with typically nearly constant 
thresholds (see e.g. [20|, but with increasing, decreasing 
or random transition thresholds (see e.g. [52| for titin 
macromolecule and [321 fo'^ artificial elastomeric protein) . 

In the following section we discuss this issue and pro- 
pose an extension of the proposed model able to repro- 
duce the hardening effect observed 6.17. in titin macro- 
molecule unfolding. 



IV. UNFOLDING ENERGY HIERARCHY 



The observed variability in the unfolding experiments 
(force plateaux, hardening, softening) can be addressed 
to a hierarchy of unfolding energies of the crystals, due 
to inhomogeneity effects with variable properties of the 
crystal domains p, [5l|, [53] ■ Indeed the experiments show 
an inhomogeneity of unfolding with a variable bond- 
breaking barriers j58| possibly due to interfacial energy 



effects [61[. Moreover another important effect can be 



due to the different orientation of the crystals in the 
macromolecule. Indeed the crystals show a directional 
deformation response of the folded domains, that, ac- 
cording with the loading direction follow different paths 
in the energy landscape, leading to different unfolding 
forces and energies [ij- Another effect, inducing hard- 
ening, is the so called n-effect (see [1J|) that, based on 
statistical considerations, addresses the observed harden- 
ing to a progressively reduced number of folded crystals 
available for unfolding in the macromolecule for growing 
elongations. 

To take care of these experimental effect in the fol- 
lowing we consider the possibility of variable unfolding 
energy of the hard domains. Thus, we first observe that, 
following the analysis of previous sections, based on the 
experimental force-displacement unfolding diagrams, we 
may estimate the fracture energy of each unfolding event 
using the relation 



^eiVu, riu) - ^e{'nu,nu + 1) = 0(n„) 



(15) 



where Q{nu) represents the variable fracture energy of 
the n„-th /3-sheet and iju ~ Lu(nu)/Lc(nu) is the strain 
corresponding to the unfolding threshold of the n„ con- 
figuration. Based on this relation we analyzed the ex- 
perimental length-force diagram for different unfolding 
macromolccules: for titin in 52 1, [3a], and [63, an d for 
TNfnAU protein from [U and Tenascin-C from |3|. The 
results are described in FigOand interestingly show a lin- 
ear empirical law 



Q{nu) = Q + (n„ - l)Ag, 



(16) 



where Q = Q(l) represents the fracture energy of the 
weakest folded domain, that has the important role of 
regulating the stability of the initial unfolded configura- 
tion and the initial unfolding length L„(l), whereas AQ 
is a fixed energy increment for successive unfolding events 
(see Fig. O. 

It is important to remark that in the case of increasing 
unfolding energies, the di-block approximation can fail, 
with the order of unfolding that is no more regulated by 
the interfacial energy effects discussed in Sect|lll In this 
case an exact solution of the problem requires a numer- 
ical analysis. To keep an analytical treatment also in 
the case of variable unfolding energies, we here suppose 
that both the mixing entropy contribution and the inter- 
facial energy term are negligible as compared with the 
unfolding energy increment AQ. Under this simplifying 




FIG. 5: Unfolding energies as a function of n„ deduced from 
the following experiments: a) AFM experiment on Titin from 
[5^ 1: b) AFM experiment on TNfnAU protein from [SI; c) 
AFM experiment on Titin from [3g]; d) AFM experiment on 
Tenascin-C from [20]; e) AFM experiment on Titin from [62 |. 



assumption, we may first easily extend the considerations 
described in Fig. [3|to obtain that the domains unfold one 
at a time in the order of their unfolding energies. Then 
we may again explicitly evaluate the unfolding lengths 
([TT[) and forces ([r5[) by simply using ([TT[) with a variable 
parameter 



£.{nxt 



8Lp 
knTl, 



■Q{nu) 



(17) 



measuring the variable ratio of dissipated and elastic en- 
ergy of the /3-sheets. 



V. AN EXPLICIT EXAMPLE: TITIN 
UNFOLDING 

To show the feasibility of the proposed model in quanti- 
tatively predicting the experimental behavior of macro- 
molecule unfolding, in this section we analyze the dif- 
fusely studied AFM stretching experiments of titin, the 
protein responsible of the passive strength of muscles. 
These proteins are very lon g ni acromolecules with con- 
tour length larger than 1 /im [60[ , whose secondary struc- 
ture is characterized by the presence of immunoglobulin 
(Ig) and fibronectin type III (FNIII) domains, folded in 
forms of /3-sheets, connected to the PEVK domain (rich 
in proline, glutamate, valine and lysine, see e.g. [30|). At 
low forces the elasticity is regulated by the tertiary struc- 
ture and the (random coil) domain orientation, combined 
with the elasticity of PEVK domains [H, Ho, [6l|. At 
higher forces the macromolecule response is dominated 
by an energetic competition of the cntropic elasticity of 
the unfolded fraction of (Ig) and (FNIII) domains and by 
the hentalpic contribution of the folded— > unfolded tran- 
sition of /3-sheets. 

As anticipated in the Introduction, one important ef- 
fect in the case of titin is the rate dependent behav- 
ior of the unfolding events (see e.g. (52|) requiring a 
MD approach as described in [20| or a reaction the- 
ory approach [18|. Here, with the aim of deducing 
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The total energy, using (flO)) . can be rephrased as 
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FIG. 6: Comparison between the AFM experiment for the 
titin protein reproduced from [S^l (continuous curves) with 
the Force vs. Elongation curves deduced by ^, pi|l . (|17|l . 
and (|16|) . Here we considered the parameters: lo = 58 nm, 
Ic = 28.43 nm, Ip = 0.36 nm, Q = 770, AQ = 420. 



a fully quantitative representative analytical approach, 
following 0, H, [2I, [55|, we take care of the observed 
rate-dependence (i.e. rate dependent unfolding energy 
barriers and dissipation), by considering effective, rate- 
dependent dissipation energies Qi. Moreover, since typ- 
ically no refolding is detected during unloading [53|, we 
assume that the (hard-soft) transitions as irreversible. 

In Fig|6]we show the ability of the model in describing 
quantitavely the behavior of titin unfolding experiments 
reported in [S^], based on the relations ^, ([TTj). and ([T7|) . 
by using the empirical law ([T5|) for the variable unfolding 
energy with the values deduced by [52l | and reported in 

FigEL- 



VI. CONTINUUM LIMIT 

In this section, aimed to a deduction of a continuum 
model for protein materials (see e.g. [l3): we analyze 
the continuum limit of the proposed model, obtained as 
a limit when n — > 00. We consider this limit at fixed 
total unfolded length that using ([3]) is given by 



(i.e. lc decreases as n grows) 
the unfolded fraction 



To this end we introduce 



n 



considered here as a (damage) continuum internal vari- 
able, with I'u G (0, 1) and i/„ = in the virgin state and 
!/„ = 1 in the fully unfolded state. The total contour 
length is then in the continuum case a function of the 
unfolding fraction that, using (|31), is given by 

Lc ~ Lc{Vu) ~ Lq+ Vu{Lc — Lq). 

Thus the damage variable v measures the change of con- 
tour length with in particular Lc = Lq in the virgin con- 
figuration (mu = 0) and Lc = Lc in the damage saturation 
of the fully unfolded state (j/„ = 1). 



knT n^ 

$tot = '^toti'^) - -^-^Lciiyu) + Qvu. (18) 
4ip 1 - 77 

where we introduced the rescaled expression 

Q ■= nQ 

ensuring that the unfolding energy Q decreases with 
growing n and so with decreasing lc- Here the deforma- 
tion variable depends on the continuum damage variable 
i^u according with the following relation 



?7 = ri{L, Vu) 



L 



Lc{Vu)' 



(19) 



The obtained framework can be inscribed in the clas- 
sical variational approach for damage known as pseudoe- 
lasticity J16l | , requiring the minimization of a damage de- 
pendent energy. We refer the reader to [Il| for a detailed 
discussion of this approach. Based on our irreversibil- 
ity assumption (no refolding) , if we indicate by L^ax the 
maximum attained assigned length we have the following 
behavior. To determine the global minimum of the en- 
ergy (jlSp . during loading (L = L„iax) we minimize both 
with respect to L and Vu ■ Minimization with respect to L 

(that is — '°3L ''^" = 0) delivers the equilibrium force as 
in ^ with the deformation variable defined in p^ . The 

minimization with respect to v (that is — g^ ''^"' = 0) 
delivers the damage as a function of the assigned length 



i^„ = i/(i) = 



1 + ^ljL-Lo 
Lc ~ Lp, 



(20) 



Interestingly in this limit we obtain a constant unfolding 
force (plateau) 



KbT U . K 




Using ((20)) we obtain that the unfolding begins (z^„ = 0) 
at 



L = Lt = 



V^Lo 



(Ve + x/2) 
whereas the fully unfolded state (i^^ = 1) is attained at 

VCLc 



L = LI 



(V?+V2)' 



The obtained unfolding behavior is shown in FigjT] (path 
OABC). 

During unloading {L < L„iax), since we neglect refold- 
ing, the behavior is again given by (|9]) with fixed damage 
that by (j20p is given by v^ = i^u{Lmax)- Different unolad- 
ing paths are shown in FigjT] {e.g. path DO is attained 
for an unloading at Vu — 0.25). 



lOOhF(pN) 




FIG. 7: Unfolding behavior of tiie continuum limit model. 
Bold lines represent the loading path, whereas continuous 
lines the unloading paths at different values of the unfolded 
fraction i/u ■ Here we used the same parameters of Fig. [l] 



important phenomenon of stretch induced unfolding of 
macromoleeules. By considering a di-block approxima- 
tion and by neglecting rate-dependent effects (possibly 
considering rate-dependent effective unfolding energies) 
we deduced a fully analytical model delivering the un- 
folding forces and lengths for the different equilibrium 
branches, all depending on the deduced main dimensional 
parameter ^. The results have also been extended to the 
case of variable unfolding energies, based on the empiri- 
cal law ([TB]) that we deduced by the experiments. Despite 
the adopted simplifying hypothesis, the deduced analyt- 
ical model shows a good qualitative (stick-slip unfolding 
evolution with regular spacing of the localized unfold- 
ing events) and quantitative agreement (see FiglH]) with 
the experimental behavior. We also showed that our ap- 
proach is amenable of the thermodynamic, continuum, 
limit useful for the extension of the model to the analysis 
of network of modular macromoleeules [4^ that will be 
the subject of our future studies. 



VII. CONCLUSIONS 

In this paper we propose, based on the time scale sep- 
aration dD), an energetic model for the description of the 
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